Transport and magnetization dynamics 
in a superconductor/single-molecule magnet/superconductor junction 
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We study dc-transport and magnetization dynamics in a junction of arbitrary transparency con- 
sisting of two spin-singlet superconducting leads connected via a single classical spin precessing at 
the frequency Q. The presence of the spin in the junction provides different transmission ampli- 
tudes for spin-up and spin-down quasiparticles as well as a time-dependent spin-flip transmission 
term. For a phase biased junction, we show that a steady-state superconducting charge current flows 
through the junction and that an out-of-equilibrium circularly polarized spin current, of frequency 
Q, is emitted in the leads. Detailed understanding of the charge and spin currents is obtained in 
the entire parameter range. In the adiabatic regime, HQ <C 2A where A is the superconducting 
gap, and for high transparencies of the junction, a strong suppression of the current takes place 
around ip « due to an abrupt change in the occupation of the Andreev bound-states. At higher 
values of the phase and/or precession frequency, extended (quasi-particle like) states compete with 
the bound-states in order to carry the current. Well below the superconducting transition, these 
results are shown to be weakly affected by the back-action of the spin current on the dynamics of the 
precessing spin. Indeed, we show that the Gilbert damping due to the quasi-particle spin current is 
strongly suppressed at low-temperatures, which goes along with a shift of the precession frequency 
due to the condensate. The results obtained may be of interest for on-going experiments in the field 
of molecular spintronics. 



I. INTRODUCTION 



Spintronics exploits the fact that an electron current 
consists of spinful carriers, with information stored in 
their spin state, which interact in a controlled way with 
their magnetic environment. Research in this field was 
pioneered in the 70s with the experiments of Tedrow and 
Meserveyi on ferromagnet/superconductor tunnel junc- 
tions as well as the experiments of Julliere^ on mag- 
netic tunnel junctions. It fully emerged in the 80s with 
the observation of spin-polarized electron injection from 
a ferromagnetic metal to a normal metal^ followed by 
the discovery of the well known giant magnetoresistance 
(GMR) effect^. Since then, most of the studies fo- 
cused on stationary magnetic states and the control of 
the electrical current by tuning the state of the magnet 
as in, e.g. GMR, see Ref. @ for a review. The study of 
the magnetization dynamics in ferromagnet (F)/normal 
metal (N) as well as F/superconductor (S) hybrids is 
much more recent and basically involves controlling the 
state of the magnet with the help of an applied electrical 
current. The main trigger was the experimental con- 
firmation of the ideas of Slonczewski 7 and Berger— that 
an electrical current may affect the state of the magnet 
via a spin transfer torque, see Ref. [9( for a pedagogical 
introduction. Indeed, a spin-polarized current of high 
enough density injected into a ferromagnet was shown 
to reverse the magnetization or generate a steady-state 
precessing magnetization in accordance with theoretical 
predictions. Many recent experiments have confirmed 
these facts, see Ref. [l(| for a review. Conversely, a mi- 



crowave driven precessing magnetization of a ferromag- 
netic layer under ferromagnetic resonance (FMR) was 
shown to emit pure spin currents, i.e. without any as- 
sociated net charge transfer, into the adjacent normal 
metal layers, see Refs. [ll| and [12[ for the theory. Such 
a spin current was only indirectly measured as an en- 
hancement of the Gilbert damping of the magnetization 
dynamics^. Recently, a similar experiment has been per- 
formed on S/F hybrids^. Contrary to the case of N/F 
hybrids, it was shown that the Gilbert damping is re- 
duced at temperatures well below the superconducting 
transition temperature, see Ref. [l5j for the theory. 

At the same time, there is a fast growing interest in 
controlling the spin orientation of single molecules or 
even a single or a few atoms in order to perform basic 
quantum operations. Recently, the control of the spin 
orientation of a single manganese atom in a semicon- 
ductor quantum dot could be achieved with the help 
of optical techniques^. At the molecular level, the 
experimental challenge is already at the level of de- 
signing a molecular junction. Single-molecule magnets 
(SMM) were recently contacted to metallic leads allow- 
ing electron-transport measurement thr oug h the m 17 i 18 
to probe their properties, see Refs. fl9ll20j | for recent 
reviews on SMMs. The access resistance due to the 
normal-contacts may however be a source of limitation 
which motivated the design of superconducting molecu- 
lar junctions. This was first done with semiconducting 
nanowires^L, carbon-nanotubes^r— and more recently 
with a single Cgo fullerene molecule^. Magnetically ac- 
tive metallofullerene molecules could also be contacted to 
superconducting leads and the proximity-induced super- 
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FIG. 1: (color online) Schematic representation of a localized 
spin (red dot) contacted to superconducting leads represent- 
ing a short S/SMM/S junction (S is for superconductor and 
SMM for single-molecule magnet), a) In the laboratory frame, 
the spin is precessing and the leads are spin-degenerate, b) 
In the rotating frame, the tilted spin is static and an effective 
magnetic field, h z — 0/2, lifts the spin-degeneracy as well as 
modifies the occupation functions of the particles in the leads. 



conductivity was studied via low-temperature transport 
measurements^. On the theoretical side, equilibrium 
properties of such Josephson junctions were extensively 
studied since a long time, see Refs. [28 33] without being 
extensive; the effect of magnetization dynamics has been 
discussed more recently, see e.g. [34 - [38j |. 

Motivated by the recent experimental breakthroughs 
of molecular spintronics we study, in the present pa- 
per, a model, first proposed in Refs. [34||35(, describ- 
ing a short (shorter than the superconducting coherence 
length) clean junction consisting of a single precessing 
spin, with precession frequency f2, connected to spin- 
singlet superconducting BCS leads, see Fig.[lji. The pre- 
cessing spin provides different transmission probabilities 
for spin-up and down quasi-particles. It is assumed to 
be the one of, e.g. a molecular magnet. The latter have 
a total spin of quite large magnitude, e.g. S = 10 for 
Fes in its ground state. For simplicity, we will neglect 
its quantum fluctuations and consider that this spin is 
classical. On the basis of such a model, our goal will be 
to compute the dc-transport properties of the junction. 
With respect to Refs. 34.35] we will consider a junc- 
tion of arbitrary transparency and, besides the dc charge 
current, we will also compute the dc spin-current across 
the junction. Even though such a set-up is still challeng- 
ing to realize experimentally the corresponding simplified 
model constitutes an interesting theoretical playground 
in order to study the combined effects of both supercon- 
ductivity and magnetization dynamics at the molecular 
scale. 

Magnetization dynamics drives the system out-of- 
equilibrium. This can be most easily seen by going to 
the rotating frame, Fig. [1J>, where the tilted molecular 
spin is static and an effective z— directed magnetic field, 
h z = fi/2, acts on the leads. Contrary to a usual mag- 
netic field, see e.g. Ref. [3!| and references therein, we 
will show in the following that h z affects the occupation 
functions of the leads translating the out-of-equilibrium 
nature of the problem. In the case of normal leads no 



dc charge current can flow through the junction. Never- 
theless, the non-equilibrium spin accumulation, i.e. the 
difference in chemical potential, Ml at T = 0, between 
spin-up and down quasi-particles, drives the emission of a 
spin current in the leads. In analogy with the case of the 
macro-magnet mentioned above, this spin-current will in 
turn damp the precession of the molecular spin eventu- 
ally leading to its complete alignment with the external 
applied magnetic field. 

Combining these out-of-equilibrium effects with super- 
conductivity in the leads brings new interesting features, 
even in the dc case 4 ^,. Indeed, integrating the junc- 
tion into a superconducting loop yields a superconduct- 
ing phase difference between the leads: (p = — 2tt&/$o, 
where $ is the total magnetic flux across the loop and 
$o = h/2e is the magnetic flux quantum. Cooper pairs 
may then be transferred across this phase biased junc- 
tion leading to the dc Josephson effect. Microscopically, 
this transfer relies on the existence of Andreev bound- 
states^. The phase dependence of the current is di- 
rectly related to the phase dependence of the spectrum 
of these subgap states. In the case of a tunnel junc- 
tion, without any magnetic impurity, the current phase 
relation (CPR) follows the well-known sinusoidal rela- 
tion: I c (<p) — Ic sin where Ic is the critical current 
depending on the normal resistance of the junction. In 
more complex cases significant departures from this sinu- 
soidal relation are known to take place, see Ref. [43| for 
a review. One well known origin for such a departure is 
the increase of the transparency, T, of the junction lead- 
ing to phase-coherent transfer of multiple Cooper pairs 
across the junction. In the ballistic limit, T = 1, and 
at low temperatures (T): I c (f) cx sm((p/2). Another 
known source of departure originates from preparing the 
system in an out-of-equilibrium state, e.g. such as by 
a microwave irradiation of the junction, see Refs. [44U45I ] 
for recent experiments and Ref. [46( for the theory. In the 
present problem, both of the these features are present. 
We will show that they strongly affect the CPR in a way 
which is proper to the (time-dependent) magnetic inter- 
face we consider. Concerning the out-of-equilibrium spin 
current we will show, in accordance with the case of the 
macro-magnet considered above, that it is strongly re- 
duced deep in the superconducting phase. This will in 
turn suppress the damping of the precession providing a 
self-consistent check of the robustness of the transport 
properties of the junction. The remaining superconduct- 
ing part of the spin current will also be shown to shift 
the precession frequency. 

The paper is organized as follows. In Sec.|TT]we present 
the mode l 34 ' 35 of a Josephson junction with a precessing 
spin. In Sec. IIIII we solve this model analytically, in the 
tunnel limit and zero temperature, to understand the ba- 
sic physics of the problem. In order to deal with a junc- 
tion of arbitrary transparency we present, in Sec. IIV1 a 
re-summation method— based on a usual Green's func- 
tion approach combined with a unitary transformation 
to the rotated frame. This method provides a convenient 
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interpretation of the out-of-equilibrium features of the 
system. It is also simple to implement numerically. In 
Secs.lVland lVII the combined analytical and numerical re- 
sults, for arbitrary transparency and finite temperatures, 
are presented and discussed for the charge and spin cur- 
rents, respectively. Finally in Sec. IVIll we summarize our 
results and conclude. 

In the following h = ks = 1 everywhere unless speci- 
fied. 



II. THE MODEL 

The model we shall consider is based on the following 
time-dependent tunnel junction Hamiltonia n 34 ' 35 : 



H(t) =H L +H R + H T {t), 



(1) 



where the perturbation theory series for the time- 
dependent perturbation, Ht, will be summed to infinite 
order in what follows. In Eq. (JXJ) , the first two terms 
correspond to the Hamiltonians of the right (R) and left 
(L) spin-singlet BCS superconducting leads: 



k,cr k 



h.c], 

where £k is the single particle spectrum, the 
(temperature-dependent) gap is given by: A a — 
A(T) e IX = and Xa is the superconducting phase in lead 
a = R,L. The last term in Eq. ([!} corresponds to the 
tunneling Hamiltonian between the leads: 

H T {t)= J2 [4,k, ff T*At) cw + h.c], (3) 
with a time-dependent transmission amplitude reading: 



T(t)=T l 2 +T 1 S(t).a, 



(4) 



where I2 is the unit 2x2 matrix and a = {<r x , <r y , <j z ) is 
the vector of Pauli spin matrices. The time-dependence 
in Eq. Q originates from the precessing motion of the 
classical spin localized in the junction. The correspond- 
ing classical equation of motion reads: 



d t S = - 7 SxH e// , 



(5) 



where 7 is the gyromagnetic ratio and H e ff = H e ffZ 
is an effective z— directed magnetic field including the 
applied field as well as other contributions such as crystal 
anisotropy and demagnetization fields. The solution to 
Eq. ^ reads: 



S(t) = S (sin 9 cos ftt, sin#sinftt, cosf 



(6) 



where 9 is the tilt angle of the spin with respect to the 
z— axis and fi = 7 -He// the precession frequency around 
this axis. 



In spin space, the transmission amplitude matrix 
reads: 



T{t) = T 1 2 



T\\a z 



Tj_a x e~ 



(7) 



where To is the direct transmission amplitude whereas 
T11 = T\S Z and Tj_ = T±S± are the spin-conserving and 
spin-flip transmission amplitudes, respectively. The lat- 
ter two depend on the magnitude, S, of the spin localized 
in the junction as well as on its orientation, 9, with re- 
spect to the quantization axis: 



T 



T s cos 9, T ± = T s sin 9, T s = T x S. (8) 



In the absence of precession, il — 0, and for normal 
metallic leads, the model in Eq. ([T]) was first proposed 
by Appelbaurn2& in order to explain tunnel conductance 
anomalies due to magnetic impurities. Its microscopic 
justification was provided by Anderson^ who has shown, 
on the basis of the Schricffcr- Wolff transformation, that: 



Ti 



2V. U 

y mix / 1 tt\ 



(9) 



where V m i X mixes the spin-states of the conduction elec- 
trons with the one of the localized spin in the junction, 
e<j < is the energy of an electron on the localized state 
and ed + U that to add a second one taking into account of 
the Coulomb charging energy U. If either or Cd + U are 
close to the Fermi level, tunneling through the magnetic 
impurity will dominate direct tunneling, i.e. T\ S> T). In 
the case of a superconducting junction, still with f2 = 0, 
the model was first considered by Kulih^S and Bulaevskii 
et ali^ for an ensemble of impurities in the dielectric 
layer between the superconductors. In the lowest order 
of perturbation theory, they have shown that, in the case 
where Ts > To, the junction is in a 7r— state. Recently, 
the model was reconsidered by Zhu and Balatsky 3 ^ for 
a single classical precessing spin in the junction. At the 
lowest order in perturbation theory, they have shown that 
the dc Josephson current is not modulated in time by the 
precession of the spin. Latter studies have focused on the 
back-action of the current on the precessing spin showing 
a possible nutation 3 ^. 

We reconsider here the model of a single classical pre- 
cessing spin in a superconducting junction keeping in 
mind that the single spin may correspond to, e.g. an 
SMM such as Mni2 or Fes- The latter have large S which 
should favor the tunneling through the spin: Ts > T\. 
Moreover, at the molecular level, one may expect that 
the tilt angle 9 may be varied in a larger range than 
for a ferromagnet under FMRiL. The cases where either 
T|| < T± (large tilt angles) or T\\ > T± (small tilt angles) 
will therefore be considered in the following. 

As stated in the Introduction, our goal is to explore the 
transport properties of such a superconducting junction: 
charge as well as spin currents, for a single conducting 
channel, arbitrary transparency and temperatures below 
the superconducting critical temperature. 
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III. THE TUNNEL LIMIT 

In this section we consider the simple limit of a tunnel 
junction. The results obtained in this limit allow for a 
direct understanding of the basic influence of the process- 
ing spin on the current flowing through the junction as 
well as the basic back-action of the (spin) current on the 
magnetization dynamics. 



A. Charge current 

In the tunnel limit, the current at lead a = R 7 L may 
formally be separated into normal an anomalous contri- 
butions: 



(10a) 



licit) 



i c a At) = c 



dt' 



dt' 



([A c a (t),A c J(t')})+h.c(jfib) 
([A c a (t),Al(t')})+ h.choc) 



where 7° G is the normal contribution to the charge cur- 
rent and J£ F the anomalous one. In the following we 
will compute the current at the left lead (a = L). The 
current at the right lead follows from charge conserva- 



tion: Ift = 



-II, so that in the rest of this section we will 



drop the lead index. In the interaction representation, 
the operator: A C L = A c , appearing in Eq. (fTU|) . is defined 



,47. 



^°(*) = E C kk, CT W T ^'WcL^'W, (11) 
k,p,CT,cr' 

and depends on the spin- and time-dependent tunnel- 
ing amplitudes which were defined in Eq. ([7]) and below. 
For a spin-singlet superconducting junction the dc charge 
current is time-independent 3 4 and, in what follows, we 
shall focus on its precession frequency dependence. 

In the absence of applied bias: Iq = 0, and only the 
anomalous part contributes to the total dc current. The 
latter reads: 

I G = 4e ( (T 2 - If) ReV R (0) - T? RcV R {Q) ) sin^, 

(12) 

where (f = \R ~ Xl, the total phase difference across the 
junction, was gauged out from the gap so that: A Q = 
A is real. The charge current of Eq. (fT"2")) depends on 
the reactive part of the anomalous 2-particle propagator 
defined (in imaginary time) as: 

v( - tuj ) = J2^J2 fe^fa te - «•>), ( 13 ) 



k,q 



from which the retarded and advanced functions are ob- 
tained by an analytic continuation: ico —> Lu±ir/ (77 = + ), 
respectively. The Matsubara superconducting Green's 



functions are defined in the usual way^ 7 -: Q a is the nor- 
mal component and T a = T\ (A € R) are the anomalous 
components of the Green's function in lead a. In the fol- 
lowing, we will mainly work with quasi-classical Green's 
functions that we simply define as: g a = ^2^G a ,k, 
fa = fa = J2kJ~a-M- For non-interacting leads they read: 



(°>(ie) 



IC 



Ve 2 + A 2 



, fW(ie) = Trv N 



A 



A 2 ' 
(14) 

where the single particle spectrum of the electrons in 
the leads has been linearized around the Fermi surface, 
vn is the normal state density of states and A = A(T) 
is temperature-dependent. With the help of Eq. (Ti"4")) . 
Eq. CES]) reads: 



1 



A 2 ] [(„• 



A 2 



where the integrand^ is a pure branch-cut. This is a 
consequence of the singular BCS density of states at the 
gap edges. These gap edge singularities will affect the 
response function, a signature of the fact that extended 
states contribute to the current. As a result, the super- 
current will have a non-analytic dependence on the pre- 
cession frequency. This can readily be seen from the zero 
temperature (T = 0) expression of the propagator (the 
reactive part being even and the dissipative part odd with 
respect to x we only consider x > in what follows): 



ReV R (x) 
lmV R (x) 



.,2 A 



AK(x), 



'N 



N 

A 



x + l 



K 



\x + lJ 



X < 1, 
X > 1, 

e(x-i) 



(15a) 



(15b) 



where x = w/2A, K is the complete elliptic integral of 
the first kind and 9 is the Heaviside function. The elliptic 
integral is logarithmically singular for ui = 2A which, by 
Kramers-Kronig, is related to the finite jump of ImT> R 
at this frequency. With the help of Eqs. (Ti"2"j) and (|15a[) 
the charge current reads: 



r = 



2eA((i 2 -i 2 )-t 2 L f sinp, fi < 2A, 



2eA l(tl-tl)-t\^K(^) sin^, Q > 2A, 



(16) 

where reduced hopping amplitudes have been introduced: 
ti = Ti/W (W = l/irv N is the band-width). 

As first noticed by Kulik^ and Bulaevskii et al^, for 
a static spin, the results of Eq. (TTB)) first show that the 
presence of the spin in the junction reduces the current. 
As will be shown by a calculation to all orders in the 
tunnel amplitude, in Sec. [V] a crossover to a ir— junction 
indeed takes place when the spin-conserving, t\\, and/or 
spin- flip, tx, tunnel amplitudes become of the order of 
the amplitude, t , for direct tunneling. 

Moreover, at zero precession frequency, the current due 
to to or tu and t± corresponds to the usual tunnel limit 
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FIG. 2: Lowest order diagram for the Josephson current 
showing the transfer of a Cooper pair from one lead to the 
other. The curvy lines correspond to the absorption/emission 
of a quantum of precession (MY) while e CT and h a denote par- 
ticles and holes of spin <r, respectively. This diagram shows 
that the Josephson current is of second order in fl/2A. 



external source, an infinite number of extended states 
are available to carry the current^. 

The regime in which Eqs. (fT8|) and (fl9|) are valid re- 
quires high precession frequencies, at least of the order 
of the amplitude of the superconducting gap. It may be 
difficult to reach, by an order of magnitude, in practice 
and would lead to a heating up of the system as wit- 
nessed by the fact that transport in this regime is domi- 
nated by extended states. In the following, we will there- 
fore be mainly interested in the adiabatic regime where: 
SI -C 2A, and the electronic degrees of freedom (of time 
scale fi/A) adjust instantaneously to the magnetic ones 
(of time scale 2tt/VL). 



of the Josephson current for a 0— or 7r— junction, respec- 
tively. The situation becomes more interesting when the 
system is driven out-of-equilibrium by the external classi- 
cal source and the spin in the junction precesses. As can 
be seen from Eq. (|16[) the precession frequency dependent 
part of the super-current is entirely carried by the spin- 
flip term. In order to single out this contribution we will 
therefore focus, in what follows, on: 8I C = I C (Q) — I c (0). 

At low frequencies, Cl < 2 A, the charge current in- 
creases with increasing precession frequency (recall that 
fl > 0). In particular, in the adiabatic limit (f2 <C 2A), 
the supercurrent arising from spin-flips reads: 



sr 



et 2 , A 



^) sin^, ft«2A, (17) 



where the lowest order term is of second order in the 
precession frequency. This comes from the two-particle 
nature of the Josephson current, see Fig. [2] Indeed, in 
order to transfer a Cooper pair from one lead to an- 
other, an Andreev spin-down/up electron together with 
its retro- reflected spin- up/down hole will flip their spins 
absorbing/emitting a quantum of precession. The spin- 
up band is therefore separated from the spin-down one 
by an energy interval equal to il. As will be shown in 
more details in Sec. IIV Al corresponds to an effec- 
tive z— directed magnetic field applied to the leads. This 
agrees with interpreting f2 as a Zeeman splitting between 
the spin-bands, see Fig. QJ>. 

On the other hand, at high frequencies, f2 3> 2A, the 
supercurrent arising from spin-flip processes decreases: 



2A 

6I C « -2et ± A — simp, 



n > 2A. 



(18) 



Finally, in the intermediate range, a resonance appears 
at — 2A where the spin-flip super-current diverges 
logarithmically: 

let 2 / fl \ 

± A togM— -1|J sinp, fi«2A. (19) 

Following the discussion at the level of the 2-particle 
propagator, such a singularity translates the fact that, 
when the system is driven at the frequency 2A by the 



B. Spin current 

The spin current at a given lead^ can be calculated 
along the same lines as the charge current (recall that 
H = 1 unless specified) : 



i s (t) = iut) + r F (t), 
-f 

2 J-oo 

I rt 

2 



dt' 

i 

dt' 



(20a) 

([A s (t),A^(t')})+ h.c](20h) 
([A s (t),A c (t')])+ h.c], (20c) 



where 1q is the normal contribution to the spin current, 
I S F the anomalous one. Contrary to the case of the charge 
current, there is no conservation law related to the spin 
of the electrons in the leads. Rather than transferring 
spins from one lead to the other, the precessing single- 
molecule magnet ejects spins in a symmetric way in the 
two leads (Ifj = 1 S L ), see Ref. [l^ where this point was 
emphasized. Therefore, there is no global spin current 
flowing through the junction. In Eq. (l2Tfl) the operator 
A c is the one which was defined in Eq. (fTTj) and, still in 
the interaction representation, the operator A s is defined 
as: 

A S W = J2 c kk, CT (*)^ CT '^', CT "(*)cL, P>CT »(*)- 

Both the normal and anomalous parts contribute to 
the total dc spin current which reads: 



I s = 



2T 
2T\\T_ 

2T T_ 



ImD R (n)cos<p-ImQ R (n) , (22a) 
lmQ R (Q) cos(fit) - (22b) 

-(ReQ R (n) - ReQ R (0)) sin(Ot)] + 
ImV R (fl)sm(nt) + 



2T||Tl 



+ (ReV R (n) - ReV R (0)) cos(fit) 
ImV R (fl) cos(Ot) - 
~(ReV R (n) - ReV R (0)) sin(fit) 



sin ip + 



cos ip, 



6 



and the y— component is derived from the x— component 
with the help of the substitution: fit — > Qt — 7r/2, in 
This general result shows that the precession of the 
spin in the junction transfers spins into the adjacent leads 
(see below for the difference between normal and super- 
conducting leads). Moreover, Eq. (|22b|) shows that, con- 
trary to the charge current, the spin current is circularly 
polarized in the xy— plane and rotates in time at the pre- 
cession frequency of the spin localized in the junction. 
Finally, the expression of the spin current involves not 
only the reactive part of the anomalous propagator but 
also its dissipative part. This is related to the fact that 
the spin current arises from breaking singlet pairs. Con- 
comitantly, the normal part also contributes to the spin 
current and is associated with Q, the normal 2-particle 
propagator, defined (in imaginary time) as: 



k,q 



(23) 



where Q is the normal component of the superconduct- 
ing Green's function and the expression of the related 
non-interacting quasiclassical function has been given in 
Eq. (|14l) . With the help of this equation, and at zero tem- 
perature, a direct computation of Eq. ([23]) yields (x > 0): 



ReQ R (x) 
lmQ R {x) 



= ReQ R (0)-ins 2 N Aq(x), 
'x-V 



2nvj f A 
l + x 



K 



+ (x — 1)E 



x + 1 
x- 1 
x + 1 



(24a) 
(24b) 



e(x-i) 



where x = w/2A, K and E are the complete elliptic in- 
tegrals of the first and second kinds, respectively, and 
is the Heaviside function. Notice that, in Eq. (I24al) . 
the combination KeQ R (x)~ KeQ R (0) eliminates an ultra- 
violet singularity in the propagator. The remaining func- 
tion q(x) is smooth and q(x) — > in the limit x — > oo 
(i.e where A — > 0), while q(x) — > 37ra; 2 /8 in the limit 
x — > 0. Because the imaginary part has a finite jump 
at x = 1, the real part diverges logarithmically at this 
point. The fact that the normal component enters the 
expression of the spin current allows us to consider both 
cases of normal and superconducting leads. 



Normal leads 



Being particularly simple, the case of normal leads 
where A = is worth examining. In this case, the T = 
spin-current arises solely from quasi-particles and reads: 




FIG. 3: Lowest order diagram for the normal spin-current 
showing the transfer of quasi-particles of spin a, denoted as 
p a - The curvy lines correspond to the emission/absorption of 
a quantum of precession. This diagram shows that the normal 
spin current is of first order in Q/2A and that it is symmetric 
between the two leads. 



I s = 

± Z 

I s = 



?t 2 



2t\\t± 



fl cos(Oi) 



2tut± 

1°= ^fisin(nt). 



(25a) 
(25b) 
(25c) 



This equation can be written in a more compact form as: 



If 1 

■Kb* 



(26) 



where Eq. ([5|), ((6|) and (JSJ have been used and the re- 
duced hopping amplitude through the spin has been in- 
troduced: ts = Ts/W (W — 1/ttv n is the band- width). 
Eqs. (I25bj) and (|25c|) explicitly show the circular polar- 
ization of the spin current in the a;?;— plane. The x— and 
y— components of the average number of spins: 



n dt. 



(27) 



where j — x, y, emitted in the leads is zero. On the 
other hand, the z— component of the spin current is time 
independent and has a linear dependence on the preces- 
sion frequency. From Eq. (|2"T|) . we see that (If) = At 2 s . 
As schematically represented on Fig. [3l this transfer of 
spin is symmetric between the two leads and we have: 
IS. = —If, which confirms the fact that there is a non- 
zero z— component of the spin current and no associated 
charge current. These arguments show that the normal 
junction is a pure spin pump where the two periodic pa- 
rameters are the projections of the in-plane spin along the 
x and y axis, see Ref. 53( in the case of a quantum dot as 
well as Ref. [l(| where similar results are reviewed in the 
case of an F/N junction. The linearity of the spin current 
with respect to the precession frequency is related to the 
single-particle nature of the tunneling quasi-particles. 



2. Superconducting leads 

The N/SMM/N case considered in the previous para- 
graph may formally correspond to the S/SMM/S case in 
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the limit where: ^> 2A — > 0, and extended (quasi- 
particle like) states carry all the current. In the opposite 
regime of low-precession frequencies, S7 <C 2 A, the dis- 
sipative part of the propagator vanishes suggesting that 
the (spin-singlet) condensate plays a dominant role. Nev- 
ertheless, we still get a non-zero spin current arising from 
the reactive part of the normal as well as anomalous prop- 
agators. With the help of Eqs. ^ and flU} the T = 
dc spin current reads: 



4 s - 
I s = 



r = 

y 



0, 



t t ± A ( n \ 2 
— i 



cos(fM) sin^ 

t n t±A/ fl \ 2 
~4~\2A~) 

sin(fii) sin tp — 



4 V2A7 

\tj_A /fly. , 
■ 1 — ) I 3 — cos ip J sin(fM), 

t t ± A / n \ 2 



4 V2A/ 



(28a) 
(28b) 

(28c) 



- JL 4-{2a) (3-co S ^jcos(Ot). (28d) 

These equations show that, in analogy with the charge 
superconducting current, the spin-current depends on the 
superconducting phase difference and has a quadratic de- 
pendence on the precession frequency at low frequencies. 
This suggests that this spin-current is of superconducting 
nature, i.e. related to the tunneling of pairs of particles. 
Its quadratic dependence on the precession frequency im- 
plies that it is smaller than the quasi-particle spin-current 
by a factor of the order of fi/A, i.e. at least one order 
of magnitude in the adiabatic regime. Actually, because 
of the singlet nature of the pairs the average number of 
spins emitted in the leads: (J?) = for j = x, y, z, as can 
be seen with the help of Eqs. (|2T[) and (|28|) . Moreover, 
its phase dependence shows that the spin-current in the 
case of a n— junction is twice larger than in the case of 
a 0— junction. In these two cases, the expression of the 
spin current slightly simplifies and its time-dependence 
is seen to be out of phase with respect to the one of the 
quasi-particle spin current of Eq. (|2"5|) by 90°. It can be 
written in a more compact form as: 



where Wo = 
respectively. 



V 4S 2 
1 and N-k 



4 7S.H e// 



2A 

2, for tp 



7 S x H e// , (29) 
and ip = 



and has to be added as a source of damping into the 
classical equations of motion of the precessing spin: 



9 t S = - 7 SxH e// +r s 



(31) 



If the external source does not compensate for this loss, 
the damping of the precession will lead to a complete 
alignment of the spin with the applied magnetic field. 

In the case of normal leads, we see from Eq. (|2l)]) that 
the spin transfer torque corresponds to a Gilbert term 
leading indeed to a damping of the precession. The cor- 
responding Gilbert constant, at T — 0, has the following 
expression: 



as 



4*| 



(32) 



and only depends on the transparency of the junction. 
On the other hand, in the case of superconducting leads 
with ip — or tt, Eq. (|2"9")l shows that the damping is sup- 
pressed and that the superconducting spin-current only 
leads to a shift of the precession frequency: 



Q' = 9, 

where, at T = 0, Cs reads: 



Cs 



r tl n 

&=M 2A C0S9 - 



(33) 



(34) 



Eq. (|34j) shows that the shift of the precession frequency 
depends on the precession frequency itself as well as on 
the tilt angle of the precessing spin and the supercon- 
ducting phase. 

In agreement with the results of the previous para- 
graph, we arrive at the conclusion that a spin precessing 
between superconducting leads is only weakly perturbed 
by its environment with respect to a spin precessing be- 
tween normal metallic leads. Uniting these two descrip- 
tions within a common framework, based on the two-fluid 
model, this dramatic reduction of the dissipation is nat- 
ural since quasi-particles are frozen deep in the super- 
conducting state and that the main current-carriers are 
spin-singlets. 



C. Gilbert damping 

In the last paragraph, we have shown that a precess- 
ing molecular spin transfers a spin current to its environ- 
ment (either metallic or superconducting leads), even if 
no voltage bias is applied. Following the Introduction, 
the emission of a spin current by the precessing spin in 
the junction implies in turn that it is loosing angular 
momentum. The corresponding spin transfer torque^ is 
given by: 

f° = II - If, = 211, (30) 



IV. GENERAL TRANSPORT EQUATIONS 

In order to go beyond the tunnel limit, we present in 
this section a unitary transformation which relates the 
time-dependent problem to a time-independent one with 
non-colinear magnetization. This will allow us to reduce 
the Dyson equations for the dressed Green's function to 
a set of algebraic equations. In particular, an explicit ex- 
pression for the Keldysh Green's function can be obtained 
from which the charge and spin currents may easily be 
computed, at least numerically. 
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A. Unitary transformation 

In order to take into account of both spin-flip processes 
and superconductivity in a compact way, we consider a 
space which is the tensor product of spin and Gorkov- 
Nambu spaces, denoted as spin£g)Nambu in what follows. 
In this 4— dimensional space, we introduce the following 
bispinor: 



^a,k — ( c Lk,f C a,k^' c ",-k,t: c a,-k. 



I) 



(35) 





( cj i 










Oi \ 


6i = 


\ o 


-Oi J 


Ti = ( 




o ) 



where the hat denotes operators acting in spin(g)Nambu 
space. For example, the spin and Nambu Pauli matrices, 
b and f , respectively, acting on the bispinor of Eq. (|35|) 
are defined as: 



(36) 



where i = 0, 1, 2, 3 and <tq — I2 while the 1, 2, 3 com- 
ponents of Oi correspond to the usual x, y, z matrices, 
respectively. 

With these conventions in hand, the action associated 
with the Hamiltonian in Eq. (TT|) and evaluated along an 
arbitrary (for the moment) C contour reads: 



S[^ 



J dt [X^«,k(*) idti>aMt) - H(t)] . (37) 



In Eq. ([57]) the Hamiltonian reads: 

H(t) = ik( ( )^,ki,k(f) + 



a, k 



k,q 

where H a k is the matrix Hamiltonian of the BCS leads: 



£k io y A a 
-io y A* a -£ k 



(38) 



and T(t) is the matrix tunnel amplitude: 

f(t) { m ° 

Y -T*(t) 

where T(t) = T Q +T ll a z +T^a x e- ta - nt , see Eq. Q. Defin- 
ing the time-dependent unitary matrix as: 



(39) 



U(t) 











(40) 



where the 4x4 Pauli spin matrix 173 was defined in 
Eq. (|3"6")l . the field operators transform as: 

^kW =^,k(*)W(*), I?a,k(t) =Z^(*)^a,k(*). 



It is then straightforward to check that this unitary trans- 
formation leaves the BCS Hamiltonian, see Eq. (|38|) . in- 
variant while the transmission matrix, see Eq. (|39j) . be- 
comes time-independent: 



f = U\t)f{t)U{t) EE f(fi = 0). 



(41) 



The unitary-transformed action therefore corresponds to 
the one of a static tilted spin with an effective magnetic 
field, h Z) acting on the leada^: 



dt 



[J2^U (idt-h z o z )^-H\. (42) 



In Eq. (|42[) the transformed Hamiltonian is time- 
independent and equals the Hamiltonian of Eq. (fTJ with 
zero precession frequency. The latter has been absorbed 
in the effective magnetic field: 



h '=2> 



(43) 



acting on the superconducting leads. From Eq. ([4"3"]). we 
recover the fact, already noticed in Sec. IIII1 that the 
precession frequency corresponds to an effective Zeeman 
splitting energy. 

Gauging out the precession frequency in the transmis- 
sion term is equivalent to the usual transformation of 
going from the laboratory frame, where the spin is pro- 
cessing, to the rotating frame where the spin appears to 
be static, see Fig. [T}d. The transformation is related to 
the peculiar harmonic dependence of the spin-flip trans- 
mission amplitude on time, see Eq. ([7]) , and implies that 
the problem has a steady-state solution. It allows us to re- 
place the time-dependent problem by a time-independent 
one at the expense of dealing with a static problem with 
non-colinear magnetization. Indeed, the effective mag- 
netic field in the leads, h z , is oriented along the z— axis 
whereas the spin in the junction is tilted in an arbitrary 
direction fixing the amplitudes Tn and T±. However, the 
system is out-of-equilibrium no matter what frame one 
considers. In particular, in the rotating frame, the out- 
of-equilibrium nature of the problem is related to the 
non-trivial action of the effective magnetic field on the 
Green's functions of the leads, i.e. with respect to what 
a usual magnetic field would do. We will show this in the 
next paragraph and end this one with a simple relation 
between the Green's functions in both frames. 

The Green's functions are matrices in spin(g)Nambu 
space and are defined with the help of: 



(44) 



where 7c is the time-ordering operator, for the Grass- 
mann fields, along the C— contour and a, ft = R, L are 
lead indices (see Eq. (|48|) with h z = for the form of 
the matrix Green's function). The unitary transforma- 
tion then allows the following correspondence between 
the Green's functions in both frames: 



(45) 



9 



B. Consequences for the leads 

Before using the unitary transformation to compute 
the transport properties of the junction we start by focus- 
ing on its consequences at the level of the leads (Green's 
functions, gap equation and thermodynamic quantities). 



1. Green's functions 

To deal with the steady-state out-of-equilibrium situa- 
tion we are considering we will use the standard Keldysh 
technique, see e.g. Refs. [55[ for reviews, and compute 
the Green's functions on a closed time or Keldysh, Ck, 
contour. This procedure leads to a doubling of the de- 
grees of freedom by introducing ip+ fields on the upper, 
forward, branch and ip_ fields on the lower, backward, 
branch. Therefore, besides the usual retarded, g R , and 
advanced, g A , Green's functions defined in Sec. there 
will be an additional component: 



(t,0 = -i(Te K ^ + , a m^(t')) 



(46) 



which is the lesser Keldysh function. In what follows, 
we will work with g~ , keeping in mind that it is related 
to the standard Keldysh Green's function, g K , with the 
help of the following identity: g K = g R — g A + 2g~ . 

The unitary transformation of Eq. (|45p will there- 
fore not only affect the spectral Green's functions (re- 
tarded and advanced) but also the Keldysh Green's func- 
tion which contains information about the occupation of 
states and hence the out-of-equilibrium properties of the 
system. For non-interacting leads the lesser component 
reads: 



"(*, = -5 a ,p n F o (g£)X - gW A ) (t, t% (47) 



where n F is the Fermi occupation function in lead a, the 
retarded and advanced functions were defined in Eq. (|14|) 
and the symbol o implies time convolution. 

In spin<X>Nambu space, the result of the unitary trans- 
formation on the spectral quasiclassical matrix Green's 
function reads: 



7 a (iu) = 



g a {iu+) 
g a {iu-) f a (iu-) 

\ -/t(iw+) 






g a (iu> + ) J 
(48) 

where iuj± — ilj ± h z . At the level of the advanced and 
retarded Green's functions, Eq. (|48p shows that the ef- 
fective magnetic field acts in a way similar to a usual 
magnetic field. 

The unusual nature of the present magnetic field ap- 
pears from the effect of the unitary transformation on 
the Keldysh component of the Green's function: 



(f: )R -^ )A )(t,n (49) 



and equivalently for the function g K . Eq. (l49l) shows 
that, contrary to what happens with a usual magnetic 
field, the occupation function is also affected by the trans- 
formation. The latter becomes a matrix in spin<8)Nambu 
space and reads: 



i(t) = n F (t)U^(t) = n F (t)e 



(50) 



Though the system is in a steady state, as it has been 
shown in Sec. lIV Al the very fact that the time-dependent 
unitary transformation modifies the occupation functions 
reflects the out-of-equilibrium nature of the problem. For 
example, we recover the fact that the spin-accumulation 
(the difference in chemical potential between spin-up and 
spin-down particles): 



de Tr 



ahp(e) 



(51) 



is non-zero only in the z— direction and, at T — 0, has an 
amplitude: /J. S Z (T = 0) = 0, corresponding to the split- 
ting between spin-up and spin-down bands, see Fig. [Tp. 



2. Density of states 

In the rotated frame the BCS density of states (DOS) 
of a given lead subject to the effective magnetic field 
reads: 



\eTh z 



(52) 



Such a DOS is similar to the one determined by Tedrow 
and Meserveyi and has been schematically represented 
on Fig. [Tp. In particular an effective (precession- 
frequency dependent) gap of 2A — 2h z between spin up 
and down quasi-particle bands appears. For Q — 2A the 
spin up and down quasi-particle peaks overlap. This is 
another manifestation of the fact that the nonanalyticity 
of the current at il = 2 A, see Eq. (|19l) . originates from 
extended states. 



3. Gap equation and thermodynamic properties 

Another proof of the unusual nature of the effective 
magnetic field is that it has no effect on the supercon- 
ducting gap A. This can be seen from the gap equation, 
A Q = Vq f a (t — 0) where Vb is the effective attractive 
coupling, which involves the equal time component of the 
anomalous part of the quasiclassical BCS Green's func- 
tion. The gap equation is therefore invariant under the 
unitary transformation. The effective magnetic field we 
are considering does not affect superconductivity in the 
leads in a way similar to a usual magnetic field. 

More generally, the thermodynamic properties of the 
leads involve energy integrations over products of the 
DOS and the occupation function. Changing variables 
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of integration therefore eliminates the effective magnetic 
field and makes these quantities invariant under the uni- 
tary transformation. 

C. Transport equations 

With the help of the unitary transformation of 
Sec. IIV A[ we derive in this paragraph general transport 
formulas valid for a junction of arbitrary transparency. 
To this end, we will follow the standard approach of com- 
putin g th e current from the Keldysh Green's function, see 
Refs. |56ll57| for some references on the subject. 

In the compact spinCgNambu space the charge and spin 
current operators are defined as: 

I c (t) ^ k (t)o- T(t)^ jq (t), 

k,q 

I S W-^E A^ s f(t)^(t) (53) 

k,q 

where a s — trfi 03) give the three components of 
the spin current, the bi (i — 0, 1, 2, 3) were defined in 
Eq. (|3"o| and erf = c>2<xo. Taking the average of the cur- 
rent yields: 

I c (t) = eTr[a f(t)g^ R (t,t)], (54a) 
I'{t) = ~T C [dff{t)gl R (t,t)], (54b) 

where g2 R (t,t') is the lesser Keldysh Green's function 
defined in Eq. (|46|) and the currents are, a priori, time- 
dependent. As in Sec. IIII1 the current is computed at a 
given electrode and the lead index is dropped for simplic- 
ity. 

The computation of the current in Eqs. (|54p therefore 
reduces to the computation of the Keldysh Green's func- 
tion. Treating the transmission amplitude T as the per- 
turbation the latter satisfies the following Dyson equa- 
tion: 

g-(t, t') = [1 + g R o f] o g(°>- o [1 + f o g A ](t, t% (55) 

where the free lesser Keldysh Green's function was given 
in Eq. (pTT|). In Eq. (1551) the check denotes the fact that all 
matrices are 8x8 matrices in the left-right<g>spin<g>Nambu 
space: 



f 9RR 


9RL \ 




i) 


V 9LR 


9th ) 







and symmetric contacts have been assumed. The dressed 
retarded and advanced Green's functions, g R ' A , also sat- 
isfy their own Dyson equation: 

g R ^ (t,t') = g^ R ^ (t, t 1 ) + gW R W f o g R W (t, t'), 

(57) 

where g^ R W are given by the proper analytic continu- 
ation of Eq. (fT3j). 



In the dc limit, that we are interested in, Eq. (|55|) may 
be further simplified by assuming that both electrodes 
have identical occupation functions: rip = rip does not 
depend on a — L, R, as no bias is applied. This yields 

g-(t,t')= -(g R on F -n F og A )(t,t') (58) 

- g R o(n F of~fon F )og A (t,t'), 

which shows that the computation of the Keldysh Green's 
function further reduces to the computation of the re- 
tarded and advanced functions. Notice that, because 
T(t) is locally time-dependent, the second term in 
Eq. (|55|) is non-zero for the present problem. 

The coupled Eqs. ([55)1 and ([57)1 are integral equations 
which arc quite complicated to solve in general. Remark- 
ably, the present problem considerably simplifies by going 
to the rotated frame. With the help of the unitary trans- 
formation of Sec. IIV Al the Dyson equation for the lesser 
Keldysh function becomes: 

g-(t,t') =-(g R oh F -fi F o§ A )(t,t') 

-g R o (fi F f - ffip^j o § A {t, 0, (59) 

where T is time-independent (and real) which eliminates 
a time-convolution. The non-commutativity between the 
transmission amplitude matrix and the occupation func- 
tion is preserved because the latter becomes a matrix, 
n F , in the rotated frame. Similarly, the Dyson equa- 
tions for the retarded and advanced functions, Eq. (|57l) . 
become: 

§ R W(t,t') = f°W A \t,t') + ~g^ R ^To'g R ( A \t,t'), 

(60) 

and another time convolution has been eliminated. This 
implies that all Green's functions in the rotated frame 
depend on the difference of their time arguments which 
allows us to go to Fourier space. Extracting the left-right 
component of Eq. (JSHJ) yields: 

fflfiM = - (ffLfl( w )^F(w) - n F {u)g A R (u)} (61) 

-9lr{u) {n F {u))f - fn F {u)) IzaM 

~3ll(u) (n F (u)f - fn F (uj)^ g RR (co), 

where n F is given by Eq. ([5U)) and T by Eqs. (|4"T1) and 
(131)1) . Notice that the last two terms in Eq. (j5Tj) are non- 
zero only because the system is out-of-equilibrium, i.e. 
the precession frequency is non-zero. Formally, they ap- 
pear as collision terms where the occupation function is 
brought away from its equilibrium value. 

Similarly, the Dyson equations for the retarded and 
advanced functions, Eq. (1571) . reduce to a set of algebraic 
equations which are straightforward to solve, yielding: 

^ ) M = 4% R(Arl ^ 03fl(A) H ! (62) 
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where the denominator reads: 



n (0)R(A) _ 



(63) 



the indices a, j3 — R, L and I4 is the unit matrix in 
4— dimensional space. In Eqs. (R)2"j) and (1531 . the non- 
interacting Green's functions depend on the effective 
magnetic field, see Eq. (|48|) . 

Substituting the expressions of Eqs. (1621) in Eq. (f6Tj) 
yields the lesser Keldysh Green's function. The latter has 
then to be substituted in the expression of the current in 
order to compute the transport properties of the junction. 
To this end we have to express the current, Eqs. (l54l) . in 
terms of the transformed Green's functions: 

/C = 6 / £ Tr t^ M M], (64a) 

IS{t) = ~\j ^[hQTglRM]- (64b) 

Eq. (I64al) shows that the charge current is actually time- 
independent. On the other hand, Eq. (|64b|) . shows that 
the spin-current is circularly polarized in the xy— plane 
whereas the z— component is time-independent. This 

comes from the fact that: <r s {t) = (5i(i), ^(t), (T3), 
where: 



Hit) 



(65) 



These results agree with the results of Sec. IIIII 

Formally, Eqs. ([Ml) together with Eqs. flEE]), (JBSJI as 
well as the self-consistent gap equation, allow for an ex- 
act computation of the charge and spin currents at any 
temperature, value of the precession frequency and trans- 
parency of the junction. Practically, the analytical com- 
putation of the current in the transparent limit is ex- 
tremely tedious in the general case. As it will be shown 
in Sec. El analytic calculations can easily be performed 
for a non-tilted/static spin and may be significantly sim- 
plified in the case where the spin has a large tilt angle 
9 — » 7r/2, see also Rcf. 37.38]. In the general case, how- 
ever, these equations have to be solved numerically. It 
turns out that, within the present reformulation and sim- 
plification of the problem, the numerical implementation 
is quite straightforward. Indeed, the system of equations 
above has the same form as the one for an equilibrium 
problem. Out-of-equilibrium effects simply appear as ef- 
fective magnetic field dependencies of the non-interacting 
Green's functions, the matrix occupation function as well 
as the sigma- matrices of the spin-current. 



V. RESULTS FOR THE CHARGE CURRENT 

The main quantity of interest in this section is the 
charge current which has been computed in the tunnel 
limit in Sec. IIIII and for which a formal exact expression 
in the transparent was given in Sec. IIV CI see Eq. (|64a[) 
together with Eqs. and (|^|) . We will start from the 



simplest cases (zero precession frequency, large tilt angle) 
before going to the general case where a fully numerical 
approach is required. Besides the current itself, we will 
also be interested in the current-carrying states, the na- 
ture of which will help us interpret the results obtained 
for the current. Information on these current-carrying 
states, as well as their occupation, is contained in the 
(charge) current kernel: 

J C M = eTr^of^H] , I C = J^ (66) 

By definition, the poles of the current kernel, or zeros 
of Eq. (|63]l (when they do exist), correspond to well- 
defined bound-states with energies below the supercon- 
ducting gap. The current due to these Andreev states 
may formally be written as: 



T c — 

1 ABS — 



+A{hz) duo 



A(h z ) 



2tt 



(67) 



where the limits of integration may be affected by the 
precession of the spin. On the other hand, the branch 
cut structure of the current kernel implies that extended 
states may carry the current as well, see Sec. IIIII where 
this fact has already been discussed in the tunnel limit. 
The current due to these extended may be formally writ- 
ten as: 



T c = 

ext 



-A(hz) , 



+ 00 



T-rW (68) 



such that the total current: I c = I% B g + ^ext- I n general, 
localized and extended states compete in order to carry 
the current^. This will clearly be seen at the level of the 
current and its dependence on the precession-frequency 
and/or superconducting phase difference. The goal of 
this section is to discuss these facts. 

In the following, numerical simulations are performed 
with a small energy relaxation rate that takes into ac- 
count phenomenologically the damping of the quasi- 
particles due to inelastic processes, such as e.g. electron- 
phonon interactions, in the leads. We will use 77 = 
10 _3 A, which agrees with typical estimates^ for usual 
superconductors, where A corresponds to the T = am- 
plitude of the superconducting order parameter. More- 
over, we will run the simulations at low temperatures, 
T = 10~ 3 A, unless specified. 



A. Static spin in the junction (Q, = 0) 

In the case where the spin in the junction is static only 
the first term in the expression of the Keldysh Green's 
function, Eq. (|6ip , gives a non-zero contribution to the 
current. Some matrix algebra yields a simple analytic 
expression for the current kernel: 



2e (To — 7s) A 2 sin<^ n^M Im 



(69) 



1 



(lj + ir)) 2 — L<J 2 (tp) 
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FIG. 4: (color online) (a) The spectrum of Andreev bound- 
states and (b) the charge current for a static spin with to = 
1 — ts and ts varies from (thick solid black line) to 1 (thick 
solid red line) with a step of 0.1. 



where r) = + , the transparencies of the junctions are 
defined as: 



To 



U 2 



(i + *§-*!) 2 + 



4t|' 



4*| 



(l + tg-i|)2+4t|' 

and, due to the spin-rotational invariance of the static 
problem, ts combines both spin-conserving as well as 
spin-flip tunneling amplitudes: t| = t 2 + 4- From 
Eq. (l69|) . we see that the current kernel has two poles 
at uj±(tp) — ±w(ip) where: 

w(p) = A y/l - % sin 2 (ip/2) - Ts cos 2 (^/2). (71) 

These poles correspond to bound-states with energies be- 
low the superconducting gap, uj{(p) < A, see Fig.[4ji. The 
integration of the current kernel reduces to an integration 
over these poles which therefore carry all the charge cur- 
rent across the junction. As can be seen from Eq. (IFj^l) . 
the Andreev bound-states carry current in opposite di- 
rections: 



4 



dip UF !■ ^ ' 



(72) 



where, because of the absence of any current-carrying 
extended state, the total current is given by: I c = Pr + 
/£.. The Fermi function then guarantees that, at low 
enough temperatures, only the lowest sub-gap state is 
occupied yielding a non-zero supercurrent: 
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FIG. 5: (color online) Density-of-states (DOS) of spin- up par- 
ticles in the rotated frame and charge current kernel (CCK) 
for: to = 0, ts = 1 and h z = 0.25A. The different line styles 
(and colors) refer to different values of the superconducting 
phase difference, tp, see the insets. The four panels correspond 
to 9 = ir/2 for (a) and (b) and = 7r/4 for (c) and (d). 



minimum in energy at (p = 0, which therefore corre- 
sponds to its ground-state phase difference. The cur- 
rent of such a 0— junction may be straightforwardly ob- 
tained from Eq. ([75]) and it is plotted on the thick black 
curve of Fig. [4]d. On the other hand, for to = and 
ts arbitrary, the spectrum of the Andreev bound-state: 

= 44/(1 +t 2 ) 2 , 
7r, see the thick 



w(» = A ^1-75 cos 2 {ip/2) with Ts 
has a ground-state phase difference of ip 
red curve in Fig. |4^,. The current of such a 7r— junction 
is the opposite of the one of a 0— junction, see Eq. ([75)1 
and the thick red curve of Fig. As was anticipated 
in Sec. IIII Al when to — ts, a transition from a 0— 
to a 7T— junction takes place. This is shown on Fig. 0] 
where it can be seen that, as soon as ts is non-zero, 
the bound-state spectrum disconnects from the contin- 
uum at ±A. Increasing ts, the curvature of the Andreev 
levels changes sign once ts > to which corresponds to 
exchanging the minimum and maximum bound-state en- 
ergies and is equivalent to a ir— shift of the phase. 



e {% - Ts) Asinv? 



tanh 



2 Jl - To sin' {ip/2) - Ts cos 2 {ip/2) 



2T 



(73) 

which is plotted on Fig. 2}). 

With these formulas in hand, we now consider the well 
known case^ where ts — and to is arbitrary. In this 
case, the spectrum of the Andreev bound-state reads: 

Lu{tp) = A yd - % sin 2 {ip/2) with To = 4*§/(l + i 2 ,) 2 . 
Assuming that we are at zero temperature, only the low- 
est in energy of these Andreev states is occupied. We 
see from the thick black curve in Fig. 0k that it has a 



B. Precessing spin in the junction {8 ~ 7r/2) 

When the spin in the junction is precessing, significant 
deviations from the results of Sec. IV Al are obtained. In 
this section, we consider the limit where the tilt angle 
of the spin with respect to the precession axis is large: 
« 7r/2. From Eq. (J5J this implies that t\\ « 0. For 
simplicity, we further assume that to — 0. Tunneling is 
then only possible via spin-flip processes through an in- 
plane precessing spin. This case has been reported on in 
Refs. |37H38j . We will summarize below the basic results 
related to the nature of the current-carrying states and 
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compare them with the results of Sec. IV Al spectrum of the bound-states. The latter correspond to 

In the case of a spin preccssmg in the plane, the sim- the poles of the current kernel of Eq. (|66l) (or equivalently 
plest quantity which may be extracted analytically is the to the zeros of Eq. (j63|). Some algebra yields: 



w(tp) 



hi 



A 2 
1-f 2 



l + f 2 cos(^-^/f 2 (l + cos^) 2 + 4(l-f 2 ) ' 



(74) 



where the transparency of the junction was redefined as: 
f = 2t\/(l + t\). Eq. (J71J) considerably differs from 
the corresponding expression in the equilibrium case, 
Eq. (T7T1) . Indeed, Eq. (1741 shows that the Andreev spec- 
trum depends now on the precession frequency via h z . 
Moreover, the appearance of the square roots suggests 
that bound-states and extended-states are inter-related. 
This is more easily seen from the asymptotic expressions 
of the bound-state spectrum: 



<<P) = 



a 11-/1,1, r^o, 



Altan^lJcos 2 ^)-^ 2 , T->1, 



(75) 

where h z = h z /A is the reduced effective magnetic field. 
Eq. ([75]) shows that, in the tunnel regime where T — » 0, 
the Andreev bound-states are dispersionless and cross 
at h z = A. This can be contrasted with the equilib- 
rium behaviour (h z = 0) where they touch the con- 
tinuum of states (this can also be seen from Eq. (fflj) 
with 7o, 71 — > 0). In both cases, however, extended 
states carry all the current (in the out-of-equilibrium case 
the currents of both Andreev bound-states cancel each- 
other). This agrees with the results of Sec. IIII Al On 
the other hand, in the transparent limit where T — > 1, 
the Andreev spectrum depends on the superconducting 
phase. As can be seen from Eq. (|75|) . the bound-states 
merge with the continuum of states for high enough val- 
ues of h z and/or a phase difference close to 7r. As has 
been shown in Refs. [37,38], this merging takes place at 
a phase difference tp c such that: duj{tp c ) / dp — 0. In 
the transparent limit, this leads to the following rela- 
tion between the critical phase and effective magnetic 
field: h Z)C {ip) = Acos 2 (c^/2) or equivalently tp c {h z ) — 
2 arccos( \fhz~f~K). Hence, bound-states exist for phase 
differences smaller than tp c and larger than 2-7T — ip c . For 
other phase differences the extended states carry all the 
current. Notice that, from these arguments, the average 
critical field at which bound-states and extended states 
merge in the transparent limit is given by: 



(76) 



In order to illustrate these facts on a concrete example, 
the density of states in the rotated frame has been plotted 



on Fig. [5^, for the special case where T = 1 and h z = 
0.25A. From the previous arguments, the critical phase 
associated with this value of the effective magnetic field 
is given by: tp c &s 0.677T. This is confirmed by Fig. [5^ 
which shows that subgap states exist for ip < tp c and 
that they tend to merge with the continuum of states, 
located at ±(A — h z ), as tp. approaches ip c . For tp > (p c 
(not represented on the figure) there are no more bound- 
states. The corresponding low— T current kernel has been 
plotted on Fig. [5Jd. In accordance with this figure, the 
total charge current may be decomposed in the following 
way (h z < A): 




-A+h z 



-A-h, 



A+h z 



A-h, 



2^ 



/», (77) 



where the first integral is over poles, u>(tp), corresponding 
to current-carrying bound-states whereas the last two in- 
tegrals are over branch cuts, of width 2h z around ±A, 
corresponding to current-carrying extended states. The 
integral has to be computed numerically which will be 
done in the next paragraph. However, we already see 
from Fig. 03 that, for tp sa (black curve), both An- 
dreev levels are occupied while for intermediate values of 
the phase, tp ~ 0.257T (blue curve), only a single Andreev 
state is occupied and that at larger values of the phase ex- 
tended states become the main current carrying states. 
These sharp changes in the occupancy of the Andreev 
states as a function of the phase difference will be seen as 
a strong suppression of the current in the current-phase 
relation (CPR) around ip rs in the next paragraph. 



C. Precessing spin in the junction (arbitrary 9) 



For arbitrary tilt angles Eq. (|64a|) of Sec. lIV Cl does not 
reduce to any simple analytic form and we proceed nu- 
merically. The results for the CPR are displayed on Fig. [6] 
and [7] in the limit of a ir— junction (t$ > to = 0). We 
see from these figures that the CPR is characterized by 
sharp steps leading to a strong suppression of the current 
around tp rj and tp m 2n. From Figs. [6] and [3 we see 
that the steps appear clearly for high transparencies, low 
precession frequencies (11 <C 2A) and preferably for tilt 
angles 8 w ir/4 or larger. For small tilt angles (6 — > 0) one 
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FIG. 6: Charge CPR for fixed transparencies: to = 0, ts = 1 
and different reduced effective magnetic fields, h z = h z /A, 
see the insets. The four panels correspond to: (a) — tt/8, 
(b) = 7r/4, (c) = 3tt/8 and (d) = tt/2. The current is in 
units of eA/h. 
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FIG. 7: Charge CPR for a given: h z = 0.25A and to = 0. 
The different line styles refer to different transparencies, ts, 
see the insets. The four panels correspond to: (a) = 7r/8, 
(b) 9 = 7r/4, (c) = 3tt/8 and (d) 8 = tt/2. The current is in 
units of eA/h. 



needs to go to higher il in order to single out the steps, 
see Fig. [6^,. Such a structure of the CPR originates from 
the existence of current-carrying bound-states which dis- 
appear, in favor of current-carrying extended states, at 
both higher values of the effective magnetic field and in- 
termediate values of the phase difference (in particular 
around (p ~ 7r). The strong suppression of the current is 
due to an abrupt change in the occupation of the lower 
and upper Andreev-levels the currents of which cancel 



FIG. 8: (color online) Charge CPR for h z = 0.25A. The 
different line styles refer to different tilt angles, see the insets. 
The four panels correspond to: (a) ts = 1, to = 0, T = 
1CT 3 A; (b) ts = 1, to = and T = 1(T 2 A (black) or T = 
10 _1 A (red); (c) ts = 1, T = 1(T 3 A and to = 0.25 (black) or 
to = 0.5 (red) or t = 0.75 (blue); (d) to = 1, T = 10 _3 A and 
ts = 0.10 (black) or ts = 0.25 (purple) or ts = 0.5 (red) or 
ts = 0.75 (green) ts = 0.9 (blue). The current is in units of 
eA/h. 



each-other when they are both occupied. 

In the case where 8 tt/2 this agrees with the results 
of the last paragraph and Figs.[5]a) and b) (see discussion 
below Eq. ((74|) ). This case, which allowed for some ana- 
lytic estimates, is however rather unphysical. In the more 
realistic case where 8 < tt/2, we see from Fig.[5Ji that the 
number of current steps doubles with respect to the case 
where 8 = tt/2. In the rotating frame, this feature is 
due to the Zeeman splitting of the Andreev states by the 
effective magnetic field, which is possible only for inter- 
mediate tilt angles. This was shown on Fig. [5J: for the 
special case 8 = tt/4. A doubling of the number of bound- 
states can be seen on this figure with respect to Fig. Oi 
for 8 = tt/2. The corresponding current kernel is plotted 
on Fig.[5ji. The latter shows that for ip s» (black curve) 
all bound-states are occupied. Because they carry cur- 
rent in opposite direction, the total current is zero. Upon 
increasing ip, one bound-state is emptied. This leads to a 
first sharp increase, in absolute value, of the current, e.g. 
tp R3 0.25tt (blue curve) on Fig. [5Ji. Increasing the value 
of the phase, another bound-state is emptied which leads 
to a second step in the CPR. Upon further increasing 
tp bound-states merge with the continuum and the cur- 
rent is again reduced. This doubling of the bound-state 
is special to the rotating frame but allows a convenient 
interpretation of the results for the current (which does 
not depend on the frame one considers). 

Notice that all these effects are washed out by temper- 
ature, see Fig. [5Jd where the parameters have the same 
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FIG. 10: Spin CPR as a function of the superconducting 
phase for fixed for ts = 1 and to = 0. The different line 
styles correspond to different tilt angles, see the insets, a) 
h z — 0.25A and b) h z = 0.1A. The spin current is in units of 
Ah/2. 



VI. RESULTS FOR THE SPIN-CURRENT 



FIG. 9: Absolute value of the charge critical current, i.e. 
maximum charge current as a function of the phase-difference 
for a given value of h z . The different line styles refer to dif- 
ferent transparencies, ts, see the insets (to = 0). The four 
panels correspond to: (a) 9 = 7r/8, (b) 9 — tt/4, (c) 9 = 3tv/8 
and (d) 9 — n/2. The current is in units of eA/h. 



values as for Fig. [5^, except for a raise in temperature. 
Indeed, by exciting particles from the lower to the up- 
per bound-state temperature broadens the levels and the 
sharp effects seen at very low temperatures disappear. 
The effects are also reduced in magnitude when direct 
tunneling is increased and the tt— to 0— junction transi- 
tion is approached. This is seen from Fig. [5J: which is 
valid for a n— junction (i<j = 1 > to) but where to in- 
creases from 0.25 (black curves) to 0.75 (blue curves). In 
the limit of a 0— junction (ts < to = 1), still with a non- 
zero tunnel amplitude through the SMM, Fig. |8ji shows 
that steps now appear around tp ~ tt but their magni- 
tude is considerably reduced with respect to the case of 
a 7T— junction. From the point of view of CPR the low— T 
tt— junction is the most interesting, as expected. 

In Fig. El the critical current, i.e. the maximum current 
as a function of the phase-difference for a given value of 
h z , is plotted as a function of h z . In the tunnel limit 
the logarithmic singularity of Eq. (fT9|) is seen at h z = 
1. When the transparency of the junction increases this 
singularity is broadened and shifted to lower frequencies. 
This is quite clear from, e.g. Fig.[9ji (limit where 8 — s- it/2 
and ts — > 1) where, upon increasing the transparency, 
the 2A— resonance shifts to a kink at the value of the 
average critical field of Eq. ([76]) in accordance with the 
discussion of the last paragraph. When the tilt angle is 
reduced the corresponding kink shifts to lower values of 
the precession frequency, see Figs. [3^, and[SJ:. 



The main quantity of interest in this section is the 
spin current which was computed in the tunnel limit in 
Sec. IIIII and for which a formal exact expression in the 
transparent limit was given in Sec. IIV CI Eq. (|64bl) to- 
gether with Eqs. (j?JTj) . (|62l) as well as the self-consistent 
gap equation. As was shown in Sec. IIIII and contrary 
to the case of the charge current, the spin current can 
be emitted in the leads only when the spin is precessing 
so it is a purely out-of-equilibrium effect. It is circularly 
polarized in the xy— plane at the precession frequency, 
il. Moreover, as the results of Sec. [111] suggest, the spin 
current is generally emitted for an inclined spin. A com- 
putation at an arbitrary tilt angle is therefore required. 
Nevertheless, the discussion on current-carrying states 
does not have to be repeated because the spin current 

kernel: I s (t,uj) = eTv\_a s (t)Tg2 R (io)\, has similar poles 
and branch-cuts as the charge current kernel of Eq. (|66l) . 
The spin current will therefore be carried by extended 
states as well as by bound-states: I s = I S ABS + ILtJ 
which are the same as the ones carrying the charge cur- 
rent. This is confirmed by Fig. [10] where it is seen that 
the amplitude of the spin CPR also has steps as a func- 
tion of the superconducting phase difference with strong 
suppression of the spin current around tp in the case 
of a 7T— junction. For identical values of the parameters 
the locations of these steps are the same as for the charge 
current. Similarly to the case of the charge current, sharp 
suppressions of the spin current therefore originate from 
an abrupt change in the occupation of the Andreev levels. 

The importance of the spin-current is related to its 
back-action on the motion of the precessing spin in the 
junction. For simplicity, we will fix the superconducting 
phase difference to tp = or tp — tt, which correspond 
to the ground-state phase difference of a 0— junction 
(where ts < to) or a it— junction (where ts > to), re- 
spectively. The fact that the spin-current may be car- 
ried by either extended or bound-states allows us then 
to distinguish between two different sources for the back- 
action: one from a normal fluid of quasi-particles (the 
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FIG. 11: (color online) (a) Case of normal leads (limit 
A — > 0): absolute value of the quasi-particle spin-current 
(black) and the dimensionless Gilbert damping constant (red) 
as a function of h z — h z /A. (b) Case of superconduct- 
ing leads: absolute value of the superconducting spin-current 
(black) and the dimensionless shift of the precession frequency 
(red) as a function of h z = h z /A. The different line styles re- 
fer to different transparencies, ts, see the insets. All figures 
are plotted for fixed: t = 0, 9 = vr/4, ip = ir and T = 10" 3 A. 
The spin current is in units of Ah/2. 

extended states) and the other one from the condensate 
(the Andreev bound-states) . Accordingly, the spin trans- 
fer torque will be separated into two contributions: 

z?s -?s i -?s 

T — T ext + T ABS 

r! xt = 2Il xt = asit f^ T) Sxd t S, (78a) 

rlss = m BS = CS( ^° ,7>) 7 S x H e// (78b) 

where the last equation holds for ip = or ip = 7r. 
Eqs. (|78a|) and (I78b|) generalize equations (|29|) and 
(|26|) of the T = tunnel limit, respectively. In the 
tunnel limit, the dimensionless coefficients as(ts,to,T) 
and £s{ts,to,T, tp) reduce to the expressions found in 
Eqs. (|32| and ([34]) . respectively. In the general case, 
they depend on both the transparency of the junction 
as well as the temperature. Notice that £s also depends 
on the tilt angle 8, see the cosine term in Eq. In 
the following, we will fix 9 = 7r/4, for simplicity. As 
already mentioned in Sec. 1111 B| as is the Gilbert damp- 
ing constant due to the normal fluid whereas (s is the 
precession- frequency shift due to the condensate. 

For normal metallic leads, the spin current and the as- 
sociated Gilbert constant are plotted at very low-T as 
a function of the effective magnetic field on Fig. [TTk for 
different transparencies ranging from the tunnel to the 
transparent limit. This figure shows that the spin current 
varies linearly as a function of the precession frequency 
in accordance with the results of the tunnel limit. This 
pumping of spins from the leads yields a Gilbert con- 
stant which depends on the transparency of the junction. 
This normal contribution does not depend on the phase 
difference and is only weakly affected by direct tunnel- 
ing, e.g. t ^ 0. When the leads are superconducting, 
this low-temperature (T) normal contribution vanishes 
altogether. The remaining low-T superconducting spin- 
current is plotted on Fig. Illb for a 7r— junction (ip = n 




T/T c T/T c 

FIG. 12: (color online) Case of a n— junction (ts > to = 
and ip = 7r): absolute value of the Gilbert damping constant, 
as, and the shift of the precession frequency, (s, as a function 
of temperature T for 9 — 7r/4. The different line styles refer 
to different transparencies, ts, see the insets. The four panels 
correspond to: (a) h z = 0.25 A, (b) zoom on (a), (c) h z = 
0.1 A and (d) zoom on (c). 



and to = 0) as a function of the effective magnetic field 
and for different values of ts- In the adiabatic regime 
(h z <C A), and actually up to h z « A, the supercon- 
ducting spin-current varies quadratically with h z . In ac- 
cordance with the tunnel limit results, the corresponding 
£s is then seen to vary linearly with h z . From the de- 
pendence of the spin-current on the precession frequency, 
we therefore recover the single-particle (linear in h z ) ver- 
sus two-particle (quadratic in h z ) nature of the normal 
versus anomalous currents, respectively. 

Still in the case of a 7r— junction, the Gilbert and 
precession-shift constants are plotted on Figs. [TU as a 
function of temperature, for different values of the trans- 
parency ts and two values of the effective magnetic field: 
h z = 0.25A (Fig. [Hi and [Hp) and h z = 0.1A (Fig. [H; 
and WM). We see clearly from these two figures that 
the Gilbert damping is strongly suppressed below tem- 
peratures of the order T « 0.1T C for h z — 0.25 A and 
T sa 0.075T C for h z = 0.1A, where T c is the supercon- 
ducting critical temperature. As shown by Figs. 112b and 
[12H . which zoom on the low-T region of Figs.[T2b and[l2b. 
respectively, the reduction of the damping upon decreas- 
ing the temperature goes along with an increasing shift 
of the precession frequency due to the superconducting 
spin-current. This shift even becomes the dominant effect 
at low-enough temperatures: T < 0.1T C for h z — 0.25A 
and T < 0.075T C for h z = 0.1 A. 

For completeness, a similar plot for a 0— junction is 
shown on Figs. [T3] In these figures h z = 0.25 but 
ts < to = 1 in Figs.fT^k andlTSb. whereas ts < t = 0.75 
in Figs. [T5b and[T5tL Provided that tunneling across the 
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FIG. 13: (color online) Case of a 0— junction (ts < to and 
ip = 0): absolute value of the Gilbert damping constant, as, 
and the shift of the precession frequency, £s, as a function 
of temperature T for fixed 6 = 7r/4 and h z = 0.25A. The 
different line styles refer to different transparencies, ts, see 
the insets. The four panels correspond to: (a) ts < to = 1, 
(b) zoom on (a), (c) ts < to = 0.75 and (d) zoom on (c). 



impurity still takes place, i.e. ts ^ 0, these figures show 
the same tendency as Figs. Q2] for the damping to be 
strongly suppressed below T s» 0.1T C in favor of a shift 
of the precession frequency due to the superconducting 
spin-current. The magnitudes of the constants are how- 
ever reduced by an order of magnitude upon reducing 
to- 

For a junction of arbitrary transparency, we there- 
fore arrive at the conclusion that, deep in the supercon- 
ducting phase, the precession of the spin is maintained 
though the frequency may be slightly lowered, in the adi- 
abatic regime, by the condensate spin-current. Other 
non-dissipative effects, such as a nutation^ of the local- 
ized spin due to the singlet-nature of the current carriers, 
may then affect the motion of the molecular precessing 
spin but they will not suppress the magnetization dy- 
namics. 



VII. CONCLUSION 

In conclusion we have presented a theory of charge and 
spin transport in a S/SMM/S junction where the single- 
molecule magnet is precessing at the frequency f2. The 
theory is based on a usual Green's function approach 
combined with a unitary transformation which replaces 
the time-dependent problem by a static one with non- 
colincar magnetization. Starting from the tunnel limit at 
zero temperature, both analytical and numerical results 
were presented, in a complementary way, for a junction of 
arbitrary transparency and at finite temperatures. Our 



goals were to study the out-of-equilibrium effects of the 
magnetization dynamics on the current flowing through 
the junction as well as the back-action of this current on 
the magnetization dynamics. 

For a phase biased junction, the unitary transforma- 
tion to the rotating frame has revealed that a steady- 
state superconducting charge current flows through the 
junction and that an out-of-equilibrium circularly polar- 
ized spin current, of frequency f2, is emitted in the leads. 
These currents and the corresponding current-phase re- 
lations (CPR) were then derived and the nature of the 
current carrying states was discussed. We have shown 
that the dynamics of the molecular spin in the junc- 
tion affects both the spectrum of the Andreev states and 
their occupation. At high precession frequencies, it even- 
tually leads to their disappearance as current-carrying 
states in favor of extended (quasi-particle like) states. 
At low precession frequencies, in particular in the ex- 
perimentally reachable adiabatic limit where f2 <C 2A, 
and for high transparencies of the junction, the abrupt 
change in the occupation of the Andreev states is seen 
as sharp steps in the CPRs with significant variations of 
the current. When the transmission amplitude through 
the magnetic impurity dominates the direct transmission 
amplitude (ts 3> to) and the junction is in a ir— state, 
the CPRs (both for charge and spin currents) show that 
the current is strongly suppressed around ip « and 
2tt. In the opposite case of a 0— junction, still with a 
non-zero transmission amplitude through the magnetic 
channel (to > ts ^ 0), a similar suppression takes place 
around ip w, though the magnitude of the effect is much 
smaller. Integrating this contact into a superconducting 
loop one may hope that such sharp features could be ob- 
served experimentally. In particular, our results for the 
charge CPR of a 7r— junction with a high transmission 
amplitude in the magnetic channel bear some similar- 
ity to experimental results on microwave irradiated SNS 
0— junctions. The CPR was recently measured^ in such 
a setup and the current was shown to be strongly sup- 
pressed near tp ss 7r upon increasing the rf power. 

We have also shown that the spin current has a weak 
back-action on the dynamics of the magnetization deep 
in the superconducting phase, i.e. at temperatures lower 
than the critical superconducting temperature by an or- 
der of magnitude or more. The reason is the strong 
suppression of quasi-particles well below T c and in the 
adiabatic regime where fl -C 2 A. This agrees with re- 
cent experiments on S/F hybrids^ under FMR where the 
Gilbert damping was shown to be much reduced contrary 
to the case of N/F hybrids^. We have also shown that 
the strong suppression of the damping, at low— T, goes 
along with a shift of the precession frequency due to the 
back-action of the remaining superconducting part of the 
spin-current on the molecular magnet. The fact that the 
dynamics of the magnetization is only weakly affected by 
the back-action of the current provides a self-consistent 
check for the validity of the above results concerning the 
CPR. 
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